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1  * 


Introduction 

It  has  been  shown  that  the  development  of  certain  types  of  cancer  can  be  hormone 
dependent.  Estrogens,  such  as  estradiol,  have  the  ability  to  bind  as  ligands  to  the  estrogen 
receptor  in  the  first  of  many  steps  which  could  result  in  the  activation  (agonistic  effect)  or  . 
repression  (antagonistic  effect)  of  genes  critical  in  the  mechanism  of  tumor  growth.  It  is  the 
object  of  this  study  to  relate  physical  and  chemical  properties  of  estrogen  derivatives  to  certain 
observed  biological  functions.  It  is  hoped  that  detailed  analysis  of  X-ray  crystallographic  data 
will  provide  important  information  to  assist  in  the  development  of  therapeutic  drugs.  My  role  is 
the  experimental  determination  of  the  electron  density  distribution  of  several  estrogens  as  part  of 
a  larger  study  to  investigate  a  wide  variety  of  estrogens. 


Task  1 .  Preliminary  studies  on  a  series  of  crystals  of  ‘A-  and  D-ring’  estrogen  derivatives. 

•  Develop  crystallization  methods  for  the  derivatives  which  are  not  yet  available  as 
charge  density  quality  single  crystals. 

-  I  had  previously  found  crystallization  methods  for  the  following  compounds: 

1 7^-estradiol  •  urea 
17ctr-estradiol  •  Vz  H2O 
17/?-estradiol  •  Vz  MeOH 
17/?-estradiol  •  2A  MeOH  •  'A  H2O 

The  best  conditions  found  to  produce  charge  density  quality  crystals  involve 
dissolving  17/?-estradiol  or  17ctr-estradiol  in  “wet”  methanol.  Allowing  the 
solvent  to  evaporate  as  slowly  as  possible  (~2-4  weeks)  will  yield  charge  density 
quality  crystals.  Unfortunately,  this  can  result  in  the  formation  of  different  crystal 
systems  (such  as  17jff-estradiol  •  Vz  MeOH  or  17/?-estradiol  •  %  MeOH  •  lA  H20). 
Selective  crystallization  was  never  achieved.  This  was  probably  due  to  the 
extremely  small  difference  in  packing  energy  of  the  different  systems.  “Wet” 
methanol  is  obtained  by  simply  allowing  methanol  to  sit  for  extended  periods  to 
absorb  atmospheric  water. 

Some  effort  was  made  in  the  past  year  to  obtain  crystals  of  diethylstilbestrol 
(DES).  This  molecule  was  not  listed  in  my  statement  of  work,  however  it  is  a 
molecule  known  to  interact  strongly  with  the  estrogen  receptor  an  exhibit 
agonistic  activity.  It  crystallizes  in  a  centrosymmetric  space  group  with  the 
molecule  sitting  on  a  special  position  resulting  in  only  half  the  molecule  being 
unique.  This  makes  it  a  good  candidate  for  charge  density  studies.  Unfortunately, 
slow  evaporation  studies  and  solvent  diffusion  studies  using  many  solvents  such 
as,  isooctance,  p-chlorophenol,  acetonitrile,  acetone,  and  several  different  types  of 
alcohols  were  not  successful.  Only  a  powdery  residue  was  observed. 
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•  Temperature  studies  on  each  derivative  to  establish  tolerances  and  the  appropriate 
temperatures  for  the  measurements. 

-  Temperature  studies  were  performed  on  the  four  crystal  systems  previously 
mentioned,  as  well  as  estrone  and  estriol.  There  were  no  problems  observed  with 
crystal  stability  in  any  of  the  systems  when  approaching  liquid  nitrogen 
temperatures.  All  of  the  systems  were  studied  both  in  house  and  at  the 
synchrotron  source  at  Argonne  National  Labs  in  order  to  obtain  data  below  liquid 
nitrogen  temperatures  (~20K)  using  liquid  helium.  We  encountered  problems 
with  crystal  stability  when  reaching  such  temperatures.  Several  attempts  were 
made  to  cool  the  crystals  by  flash  cooling  as  well  as  by  cooling  the  crystals  very 
slowly.  Neither  method  improved  the  crystal  stability  at  such  low  temperatures. 


•  Preliminary  routine  X-ray  crystal  structure  determination  on  each  derivative  to  check 
for  composition,  quality,  and  solvation. 

Preliminary  data  sets  were  collected  on  the  four  systems  crystallized.  Two  of 
these  structures  were  as  expected.  17/3-estradiol  •  Vi  MeOH  was  a  previously 
unknown  structure  and  was  electronically  submitted  to  Acta  Crystal )  The 
structure  of  the  17/3-estradiol  •  2/3  MeOH  •  lA  H20  system  was  also  new,  however 
while  chemically  interesting,  is  not  a  good  candidate  for  charge  density  research. 
The  extreme  number  of  atoms  in  the  unit  cell  makes  a  detailed  charge  density 
study  of  the  system  very  difficult.  This  system  is,  in  fact  very  important 
chemically,  because  the  three  independent  molecules  in  the  asymmetric  unit  are 
not  identical.  While  they  share  the  same  stereochemical  configuration,  they 
differ  slightly  in  physical  geometry.  This  structure  is  discussed  in  detail  in  the  1st 
publication  from  this  study.(2)  This  publication  is  attached  as  supplemental 
material  to  this  report. 


Task  2.  Electron  density  quality  data  collection  on  the  above  mentioned  estrogen  analogues. 

•  X-ray  diffraction  studies  at  liquid  nitrogen  temperatures  on  crystals  that  did  not 
qualify  for  lower  temperatures. 

-  A  complete  data  set  of  17or-estradiol  •  V2  H20  has  been  successfully  collected  and 
analyzed. 

-  A  complete  data  set  of  17/3-estradiol  •  urea  has  been  successfully  collected  and 
analyzed. 

A  complete  data  set  of  17/3-estradiol  •  V2  MeOH  been  successfully  collected  and 
analyzed. 

•  X-ray  diffraction  studies  at  liquid  helium  temperatures. 

-  Work  has  been  continued  to  try  to  find  a  better  strategy  for  cooling  crystals  down 
to  near  liquid  He  temperatures. 
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Task  3.  Interpretation  and  analysis  of  nitrogen  and  helium  temperature  charge  density  data  sets 
of  above  mentioned  estrogen  analogues. 

•  Analysis  of  the  experimental  data. 

-  During  this  study  it  was  discovered  after  attempting  multipole  refinements  of  the 
three  systems  that  our  process  for  data  treatment  was  not  consistent  enough  to 
yield  reliable  results.  The  problems  lie  in  the  integration  of  the  raw  data.  This 
required  that  we  take  a  step  back  and  reevaluate  the  application  of  the  software  we 
use. 

-  Integration  of  the  raw  data  involves  integrating  the  intensity  of  the  reflections  as 
measured  by  a  two-dimensional  CCD  detector.  Several  parameters  must  be 
defined  to  determine  exactly  how  the  software  integrates  the  reflections.  It  was 
the  combination  of  parameters  that  had  to  be  refined.  There  are  three  parameters 
which  critically  affect  the  outcome  of  the  integration,  and  they  are: 

-  Box  Size  -  area  on  2-dimensional  frame  to  be  integrated  for  each  reflection 

-  Profile  Fitting  Limits  -  threshold  for  reflections  which  are  used  to  determine 
the  profiles  applied  in  the  fitting. 

-  Simple  Sum  Perimeter  Limit  -  Determines  how  far  out  on  the  reflection 
profile  to  integrate 

It  was  found  that  different  detector  settings,  even  for  the  same  data  set,  require 
different  box  size  parameters  and  profile  fitting  limits.  The  simple  sum  perimeter 
limit  of  0.02  was  found  to  be  the  best  value  for  all  three  data  sets.  Selected 
crystal,  integration,  and  reflection  data  for  the  three  systems  are  shown  in 
Appendix  A. 

-  Upon  completion  of  the  integration,  more  data  treatment  is  required.  This 
includes  the  absorption  correction,  data  scaling,  and  statistical  outlier 
determination. 

-  The  multipole  model  requires  a  local  coordinate  system  be  set  up  for  every  atom. 
A  standard  coordinate  system  was  developed  and  applied  to  each  structure  where 
possible.  The  submitted  version  of  the  coordinate  system  paper  is  attached  as 
supplemental  material  to  this  report.(3) 

-  It  was  found  that  the  starting  values  for  the  multipole  model  greatly  influenced  the 
path  the  refinement  would  take.  It  was  determined  that  a  specific  set  of  starting 
values  should  be  applied  to  each  structure  to  ensure  consistency.  These  optimum 
values  are  shown  in  Appendix  A.  . 

-  The  multipole  refinement  as  well  as  full  topological  analysis  was  completed  for 
all  three  structures. 

17/?-estradiol  •  urea:  As  stated  earlier  in  this  study,  the  first  data  set  was  not 
usable.  The  second  data  set  collected  was  of  high  quality.  This  allowed  the 
multipole  refinement  to  be  completed. 

-  17/ft-estradiol  •  Vi  MeOH:  Despite  the  fact  that  the  crystal  system  is  PI , 
meaning  there  is  no  symmetry  equivalent  data  which  reduces  redundancy  in 
the  data,  the  multipole  refinement  was  successfully  completed  as  well  as  full 
topological  analysis. 

17«-estradiol  •  Vi  H20:  The  water  molecule  of  this  system  lies  on  a  2-fold 
axis  of  rotation,  meaning  only  half  of  the  molecule  is  unique.  This  in  itself  is 
not  a  problem  except  that  it  was  discovered  during  the  multipole  refinement 
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that  the  hydrogen  atom  was  very  slightly  disordered.  The  position  it  refined  to 
generated  a  symmetry  equivalent  hydrogen  atom  which  created  an  H  -  O  -  H 
bond  angle  of  less  than  90°.  The  thermal  parameter  of  the  hydrogen  atom  is 
twice  as  large  as  a  typical  hydrogen  atom  of  the  system.  Due  to  the  fact  the 
disorder  is  a  result  of  a  shift  on  the  order  of  0.1  Angstroms  for  the  hydrogen, 
the  effect  of  this  disorder  is  taken  up  by  the  large  thermal  parameter  of  the 
hydrogen.  Unfortunately  this  disorder  greatly  effects  the  hydroxy  atoms  that 
are  hydrogen  bound  to  the  water.  This  complicated  the  refinement 
significantly,  however  it  was  successfully  completed  as  well  as  full 
topological  analysis. 

-  The  refinements  were  completed  using  k  and  k’  parameters  of  1.40  and  1.29 
respectively.  These  values  are  critical  because  a  hydrogen’s  isotropic  thermal 
parameter  is  directly  correlated  to  its  K  and  k’  parameters.  These  values  were 
initially  set  at  Kch  =  1 .40,  1 .29  and  Koh  -  1  -30,  1 .29  because  they  gave  the  most 
appropriate  ratio  of  thermal  parameter  magnitude,  based  on  neutron  data  of  other 
small  organic  molecules,  between  the  hydrogen  and  the  atom  it  is  bound  to.  In 
the  end  the  value  of  1 .20  was  used  for  K  because  it  is  a  much  more  commonly 
accepted  value.  This  required  that  a  great  deal  of  time  be  spent  re-refining  all 
three  systems  using  the  K  value  of  1.20. 

Comparison  of  the  results  from  the  series  of  estrogen  analogues. 

-  Initial  comparisons  of  the  multipole  models  and  the  topological  analysis  have 
yielded  expected  and  unexpected  results.  It  was  expected  that  the  core  structure 
of  the  estrogen  molecules  would  remain  relatively  unchanged  from  system  to 
system  even  with  chemical  substitutions  at  the  activity-sensitive  ends.  This  was 
found  to  be  the  case. 

-  One  major  question  that  needed  to  be  answered  was  could  such  small  features  as 
oxygen  lone  pair  densities  be  determined  on  such  large  systems.  If  they  could  be 
determined,  then  you  have  to  ask;  Would  they  be  affected  by  different  hydrogen 
bonding  schemes?  Would  the  lone  pair  density  of  the  oxygen’s  on  the  aromatic 
ring  conjugate  to  the  pi  system  of  the  aromatic  ring?  The  answer  is  that  we  can 
determine  features  such  as  lone  pair  densities  on  the  oxygen’s  (see  diagrams  in 
appendices).  It  was  found  that  each  oxygen  atom  had  two  lone  pairs  in  a  rough 
sp3  type  geometry.  This  demonstrates  that  the  lone  pairs  are  surprisingly  robust  in 
that  they  do  not  significantly  change  in  different  hydrogen  bonding  schemes  and 
when  the  oxygen  is  bound  to  aromatic  systems.  This  was  somewhat  surprising, 
however  this  type  of  analysis  is  not  well  represented  in  the  literature. 

-  Although  the  influence  of  the  lattice  is  inherent  in  the  multipole  model  of  the 
estrogen  molecule,  the  electrostatic  potential  is  believed  to  provide  a  reasonable 
representation  of  a  single  molecule  isolated  from  the  crystal  lattice.  In  these 
studies,  the  electrostatic  potential  of  the  hydroxy  groups  on  the  aromatic  rings 
were  consistent  in  both  shape  and  magnitude.  Generally  covering  nearly  the  * 
entire  lone  pair  region,  the  maximum  magnitudes  are  —20  eA'1.  One  could  argue 
that  this  is  expected  given  the  consistency  in  the  lone  pair  regions  of  the  oxygen 
atoms.  However,  the  C 17  hydroxy  groups  are  also  consistent  in  their  lone  pair 
regions  and  they  exhibit  an  interesting  trend.  The  potential  changes  significantly 
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in  both  shape  and  magnitude  based  on  the  relative  position  of  the  hydroxy 
hydrogen  with  respect  to  HI 7.  17^-estradiol  •  urea  and  17a-estradiol  •  Vi  H20 
displayed  a  shift  in  the  potential  toward  the  lone  pair  region  opposite  HI 7,  as  well 
as  increased  magnitudes  of  ~  -35  and  — 45  eA'1.  H20  and  H20’  were  virtually 
opposite  HI 7  in  17/?-estradiol  •  Vi  MeOH.  This  resulted  in  a  fairly  uniform 
distribution  of  the  electrostatic  potential  with  a  maximum  value  of  — 20  eA'1. 
Plots  of  the  electrostatic  potential  are  shown  in  Appendix  A. 

-  The  first  paper  in  this  series  discusses  the  17/?-estradiol  •  urea  system  is 
completed  and  has  been  submitted  for  recommendation  to  The  Proceedings  of  the 
National  Academy  of  Sciences. 

-  The  manuscript  for  the  17/?-estradiol  •Vi  MeOH  system  is  currently  being 
prepared  for  submission  to  the  Journal  of  the  American  Chemical  Society. 

-  The  manuscript  for  the  1 7  or-estradiol  •  Vi  H20  system  is  currently  being  prepared 
for  submission  to  the  Journal  of  Physical  Chemistry. 

•  Analyze  relationship  of  charge  density  to  receptor  binding  affinity  and  the 
chemical/biological  effects  as  related  to  breast  cancer. 

-  Presumably  there  is  a  preferred  orientation  for  the  hydrogen  of  the  C17  hydroxy 
group  in  the  active  site  of  the  estrogen  receptor.  Unfortunately,  the  protein  crystal 
structures  determined  to  date  do  not  have  enough  resolution  to  determine  this.  It 
seems  possible  that  this  is  one  mechanism  by  which  the  estrogen  ligands  could 
stimulate  different  responses  from  the  estrogen  receptor.  Binding  of  different 
ligands  would  most  likely  result  in  different  orientations  of  the  C 17  hydroxy 
hydrogen,  thereby  presenting  different  electrostatic  potentials  to  the  receptor. 

This  could  help  to  account  for  the  changes  in  binding  affinity  and  activity  of  the 
different  estrogen  derivatives. 

-  Continued  effort  must  be  made  to  analyze  more  systems  because  a  sample  of 
three  structures  is  not  enough  to  confirm  the  trends  that  are  appearing  in  the 
electron  density  and  the  electrostatic  potential. 

Key  Research  Accomplishments 

•  The  core  estrogen  structure  is  very  consistent  between  derivatives. 

•  Determination  that  it  is  possible  to  locate  lone  pair  densities  of  oxygen’s  in  such  large 
systems 

•  Determination  that  the  oxygen  atoms  of  these  systems  are  decidedly  sp3  in  shape 
despite  the  fact  they  are  bound  to  aromatic  neighbors. 

•  Determination  that  the  lone  pairs  appear  to  be  unaffected  by  completely  different 
hydrogen  bonding  environments. 

■  •  Determination  that  the  electrostatic  potential  around  the  aromatic  hydroxyl  groups  are 

consistent  despite  different  hydrogen  bonding  environments. 

•  Determination  that  the  electrostatic  potential  around  the  C 17  hydroxy  groups  changes 
significantly  in  both  shape  and  magnitude  based  on  the  relative  position  of  the 
hydroxy  hydrogen  with  respect  to  H17. 
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Kumaradhas,  Damon  A.  Parrish,  Nan  Wu  and  Elizabeth  Zhurova 
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Conclusion 

The  ability  to  determine  the  electron  density  distribution  of  larger  molecular  systems  such 
as  estrogens  has  been  clearly  demonstrated.  Observation  of  the  consistency  of  the  core  structure 
of  the  estrogen  molecule  provides  interesting  chemical  and  physical  information,  as  well  as  a 
basis  for  future  estrogen  studies.  The  robustness  of  the  lone  pairs  under  different  hydrogen 
bonding  conditions  was  unexpected.  Finally,  the  trends  in  the  electrostatic  potential  of  the  C 17 
hydroxy  groups  suggest  a  possible  mechanism  by  which  different  estrogen  ligands  can  elicit 
varied  responses  from  the  receptor.  Proteins  will  react  differently  in  their  conformational 
changes  to  variations  in  the  potential  of  the  ligand.  These  actions/reactions  will  determine  the 
downstream  outcome  of  the  protein  complex.  More  studies  of  different  estrogen  derivatives,  or 
the  same  derivatives  with  different  orientations  of  the  C17  hydroxy  hydrogen,  would  help  to 
support  this  hypothesis. 


Career/Research  Plans 


As  I  stated  in  my  original  proposal,  my  career  goals  were  to  complete  my  PhD  and  to 
pursue  a  career  in  research  crystallography,  perhaps  in  the  area  of  medicinal  chemistry  and 
rational  drug  design.  Working  on  this  project  has  given  me  the  opportunity  to  do  just  that.  I 
have  received  a  three-year  Postdoctoral  Fellowship  Award  sponsored  by  the  American  Society 
for  Engineering  and  Education  (ASEE).  This  fellowship  places  me  at  the  Naval  Research 
Laboratory  in  Washington  DC  to  provide  crystallographic  support  for  the  National  Institute  on 
Drug  Abuse  (NIDA).  I  believe  I  was  awarded  this  position  based  largely  on  my  work  with  the 
estrogen  systems.  I  would  like  to  thank  the  DOD  for  their  support  of  my  work  and  for  their 
facilitation  in  furthering  my  career.  I  would  also  like  to  say  that  I  highly  recommend  Dr.  Alan 
Pinkerton  and  his  current  graduate  students  as  future  candidates  for  this  research  award.  Dr. 
Pinkerton  and  the  University  of  Toledo  provide  an  ideal  environment  for  students  to  participate 
in  an  important  area  of  research,  which  is  in  need  of  quality,  knowledgeable  people. 
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Appendix  A 


Crystal  Data 

Chemical  Formula 

C19H28N2O3 

Temperature 

100.0(1)  K 

Crystal  Dimensions 

0.35  x  0.37  x  0.40  mm 

Space  Group 

P2i2i2i 

a 

7.9022(9)  A 

b 

9.2228(10)  A 

c 

24.5890(28)  A 

Volume 

1792.06(56)  A3  1 

Z  (Crystallographic) 

4 

Integration  Parameters 

Box  Size  (°) 

Profile 
Fitting  (I/o) 

Simple  Sum 
Perimeter  Limit 

Low  Angle 

1.5  x  1.5  x  1.0 

30/10 

0.02 

Medium  Angle 

1.2  x  1.2  x  0.8 

30/10 

0.02 

High  Angle 

1.0  x  1.0  x  0.6 

20/10 

0.02 

Reflection  Statistics  (from  SORT  A  V) 

Total  Reflections 

110999 

Rejected  Outliers 

779 

Unique  Reflections 

13187 

Average  Redundancy 

8.4 

Resolution 

1.180  A’1  " 

Completeness 

98.6  % 

Ri 

3.52% 

r2 

3.98  % 

Rw 

12.84  % 

Z  (Refinement) 

1.999 

Table  A-l .  Selected  crystal,  integration,  and  reflection  data  for  17/9-estradiol  •  urea. 

» 

'  V  ■ 
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Crystal  Data 

Chemical  Formula 

C37H52O5 

Temperature 

100.0(1)  K 

Crystal  Dimensions 

0.22  x  0.26  x  0.42  mm 

Space  Group 

PI 

a 

7.2934(1)  A 

b 

9.2795(1) A 

c 

12.3918(2)  A 

a 

89.4681(5) 

_ A 

87.8569(4) 

_ z _ 

70.7659(4) 

Volume 

791.290(33)  AJ 

Z  (Crystallographic) 

2 

Int 

egration  Parameters 

Box  Size  (°) 

Profile 
Fitting  (I /&) 

Simple  Sum 
Perimeter 
Limit 

Low  Angle 

1.5  x  1.5  x  1.0 

20/20 

0.02 

Medium  Angle 

1.2  x  1.2  x  0.8 

20/20 

0.02 

High  Angle 

1.0  x  1.0  x  0.6 

10/10 

0.02 

Reflection  Statistics  (from  SORTAV) 

Total  Reflections 

86369 

Rejected  Outliers 

33 

Unique  Reflections 

29051  H 

Average  Redundancy 

3.0 

Resolution 

1.329  A'1  ~  " 

Completeness 

91.9% 

Ri 

5.77  % 

r2 

5.34  % 

Rw 

1  15.25  %  1 

Z  (Refinement) 

1.219 

Table  A-2.  Selected  crystal,  integration,  and  reflection  data  for  17/?-estradiol  •  methanol. 
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Table  A-3. 


Crystal  Data 

Chemical  Formula 

C18H25O2.5 

Temperature 

100.0(1)  K 

Crystal  Dimensions 

0.24  x  0.33  x  0.33  mm 

Space  Group 

C2 

a 

19.0240(5)  A 

b 

7.0652(2)  A 

c 

13.3496(3)  A  ' 

_ i _ 

124.0553(10) 

Volume 

1486.57(10)  A3 

Z  (Crystallographic) 

4 

Into 

igration  Parameters 

Box  Size  (°) 

Profile 
Fitting  (I/o) 

Simple  Sum 
Perimeter 
Limit 

Low  Angle 

1.2  x  1.2  x  0.8 

40/10 

0.02 

High  Angle 

1.0  x  1.0  x  0.6 

30/10 

0.02 

Reflection  Statistics  (from  SORTAV) 

Total  Reflections 

85606 

Rejected  Outliers 

65 

Unique  Reflections 

14594 

Average  Redundancy 

5.9 

Resolution 

1.319  A'1 

Completeness 

98.2% 

Ri 

3.51  % 

r2 

3.88% 

Rw 

14.02  %  1 

Z  (Refinement) 

2.010 

Selected  crystal,  integration,  and  reflection  data  for  Ha-estradiol  •  Vz  H2O. 
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Figure  A-l.  Electrostatic  potential  isosurfaces  of  the  01  hydroxy  group  of  17y9-estradiol  • 

urea.  The  blue  surface  in  each  picture  represents +1.0  eA'1.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA'1,  B  -0.10  eA'1,  C  -0. 15  eA'1,  D 
0.20  eA'1. 


Figure  A-2.  Electrostatic  potential  isosurfaces  of  the  02  hydroxy  group  of  17/^estradiol  • 
urea.  The  blue  surface  in  each  picture  represents  +1.0  eA'1 .  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA'1,  B  -0.10  eA'1,  C  -0.15  eA'1,  D 
0.20  eA'1,  E  -0.25  eA' ,  F  -0.30  eA'1,  G  -0.35  eA'1. 
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Figure  A-3.  Electrostatic  potential  isosurfaces  of  the  01  hydroxy  group  of  17/?-estradiol  • 

methanol.  The  blue  surface  in  each  picture  represents  +1.0  eA'1.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA'1,  B  -0.10  eA'1,  C  -0.15  eA'1,  D  - 
0.20  eA'1. 
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Figure  A-4.  Electrostatic  potential  isosurfaces  of  the  01’  hydroxy  group  of  17/?-estradiol  • 

methanol.  The  blue  surface  in  each  picture  represents +1.0  eA'1.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA*1,  B  -0. 10  eA*1,  C  -0. 15  eA*1,  D  - 
0.20  eA*1. 
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Figure  A-5.  Electrostatic  potential  isosurfaces  of  the  02  hydroxy  group  of  17/J-estradiol  • 

methanol.  The  blue  surface  in  each  picture  represents  +1.0  e A'1.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA*1,  B  -0. 10  eA'1,  C  -0. 1 5  eA*1,  D  - 
0.20  eA'1. 


Figure  A-6.  Electrostatic  potential  isosurfaces  of  the  02’  hydroxy  group  of  17/?-estradiol  • 

methanol.  The  blue  surface  in  each  picture  represents  +1.0  eA"1.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA"1,  B  -0. 10  eA'1,  C  -0. 15  eA'1,  D  - 
0.20  eA"1. 
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Figure  A-7.  Electrostatic  potential  isosurfaces  of  the  01  hydroxy  group  of  1  Tor-estradiol  •  14 
H2O.  The  blue  surface  in  each  picture  represents  +1.0  eA'1.  The  red  surfaces 
represent  the  following  potentials:  A-0.05  eA'l,B-0.10  eA‘VC-0.15  eA*1,  D- 
0.20  eA'1. 
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Figure  A-8.  Electrostatic  potential  isosurfaces  of  the  02  hydroxy  group  of  17#-estradiol  •  lA 

H2O.  The  blue  surface  in  each  picture  represents  +1.0  eA"1.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA"1,  B  -0. 10  eA'1,  C  -0: 15  eA"'1,  D  - 
0.20  eA"1,  E  -0.25  eA"1,  F  -0.30  eA'1,  G  -0.35  eA"1,  H  -0.40  eA"1, 1  -0.45  eA"1. 


24 


Current  information: 

Dr.  Damon  Parrish 
Research  Chemist 

Naval  Research  Laboratory,  Code  6030 
4555  Overlook  Ave. 

Washington,  DC  20375 
Phone:  (202)404-2141 
Fax:  (202)767-6874 
Email:  dparrish@ccs.nrLnavy.mil 


organic  compounds 


Acta  Crystallographies  Section  C 

Crystal  Structure 
Communications 

ISSN  0108-2701 

A  new  estra-1 ,3,5(1 0)-triene-3,1 7/?- 
diol  solvate:  estradiol-methanol- 
water  (3/2/1) 

Damon  A.  Parrish  and  A.  Alan  Pinkerton* 

Department  of  Chemistry,  University  of  Toledo,  Toledo,  Ohio  43606,  USA 
Correspondence  e-mail:  apinker@uoft02.utoledo.edu 

Received  22  October  2002 
Accepted  21  November  2002 
Online  25  January  2003 

The  title  solvate  of  the  steroid  17/J-estradiol  (E2)  with 
methanol  and  water,  Ci8H24020.67CH40  0.33H20,  is  the 
first  E2  derivative  to  contain  three  crystallographically 
independent  molecules  in  the  asymmetric  unit.  The  three 
steroid  molecules,  along  with  two  methanol  molecules  and  a 
water  molecule,  create  a  three-dimensional  hydrogen-bonded 
system.  Three-sided  columns  are  formed,  with  the  estradiol 
molecules  aligned  lengthwise  parallel  to  (101),  and  joined  by 
solvent  molecules  at  both  hydrophilic  ends.  The  three 
estradiol  molecules  differ  slightly  in  their  ring-bowing 
angles,  Le.  the  angle  between  the  mean  plane  of  the  A  ring 
and  that  of  the  BCD  ring;  this  angle  ranges  from  7.1  to  12.2°. 

Comment 

17/3-Estradiol,  is  a  member  of  the  estrogen  family  of 
hormones.  In  recent  years,  interest  in  these  molecules  has 
focused  primarily  on  understanding  their  biological  role  in 
initiating  breast  cancer.  It  is  well  known  that  their  ability  to 
form  hydrogen  bonds  in  the  active  site  of  the  estrogen 
receptor  (ER)  influences  biological  activity.  This  ability  to 
form  hydrogen  bonds  has  been  clearly  demonstrated  in  an 
array  of  crystal  structures,  especially  in  that  of  E2,  containing 
different  solvent  molecules  or  other  hydrogen-bond  acceptors. 

When  comparing  these  structures,  the  flexibility  of  the 
hydrophobic  region  of  the  molecule,  more  specifically  the  B 
ring,  becomes  very  apparent.  Several  papers  have  already 
discussed  the  ring  bowing  of  this  molecule  (Cooper  et  ah ,  1969; 
Cody  et  ah ,  1971;  Busetta  et  ah,  1976;  Duax  et  al. ,  1979).  Weise 
&  Brooks  (1994)  took  it  a  step  further,  performing  molecular- 
modeling  calculations  on  observed  and  novel  confirmations. 
Ivanov  and  co-workers  carried  out  a  related  study  which 
included  a  larger  body  of  compounds  and  more  detailed 
analysis  (Ivanov  et  al ,  1998).  They  concluded  that  there  are 
two  possible  conformations  very  close  in  energy,  which  differ 
in  the  B-ring  arrangement,  although  the  chance  of  the 
strained-geometry  binding  is  low,  as  one  quarter  of  the  binding 
energy  (—11.9  kcal  mol”1;  1  kcal  mol-1  =  4.184  kJ  mol”1)  is 


predicted  to  be  lost  (Anstead  et  al ,  1997).  The  point  initially 
postulated  by  Weise  &  Brooks  remains,  Le.  that  the  flexibility 
to  allow  bending  or  other  conformational  changes  of  the 
ligand  in  the  receptor  appears  energetically  achievable,  and 
this  could  be  an  important  property  in  determining  their 
activity. 

Molecular-dynamics  studies  on  the  ligand-binding  domain 
(LBD)  of  the  ER  demonstrate  that  the  motion  of  the  LBD 
requires  the  A  ring  to  remain  fairly  steady  while  the  CD  ring 
retains  a  higher  degree  of  freedom.  The  range  of  motion  found 
in  that  study  (Maalouf  et  ah,  1998)  agrees  very  well  with  the 
reported  crystal  structures.  It  is  well  known  that  the  hydrogen 
bonding  of  the  ligand  in  the  LBD  causes  conformational 
changes  in  the  receptor.  This  seemingly  accommodating 
motion  of  the  LBD  and  the  ligand  supports  the  idea  postu¬ 
lated  by  Weise  &  Brooks  that  the  flexibility  of  E2,  or  of  any 
other  ligand  which  enters  the  LBD,  could  also  effect  the 
activity  of  the  ligand.  A  more  rigid  or  flexible  ligand  could 
change  the  natural  motion  of  the  ER  complex,  resulting  in  a 
change  in  activation  factor  (AF-2)  activity,  and  therefore 
possibly  effecting  co-activator  recruitment.  This  process  is  not 
well  understood;  the  flexibility  of  the  ligand  is  certainly  only 
one  of  many  physical  factors  associated  in  the  activity  of  E* 
The  structure  of  the  title  solvate,  (I),  continues  the  trend  in 
observing  the  flexibility  of  this  molecule  from  a  structural  and 
hydrogen-bonding  point  of  view,  while  the  molecule  lies  in  the 
lowest  energy  conformation. 


The  crystal  of  (I)  contains  three  crystallographically  unique 
E2  molecules  (Fig.  1).  The  B  rings  of  the  E^  molecules  adopt 
the  typical  conformation  of  a  distorted  7a,8£-half-chair,  and 
this  is  responsible  for  most  of  the  structural  flexibility. 
Calculation  of  the  ring-bowing  angle,  as  defined  by  Duax  & 
Norton  (1975),  reveals  a  range  of  5.1°  for  the  three  molecules. 
Table  1  compares  the  ring-bowing"  angles  of  the  currently 


A  view  of  molecule  A  of  (I),  with  the  atom-numbering  scheme;  the  other 
two  molecules  are  similarly  labeled,  with  the  corresponding  suffix  B  or  C. 
Displacement  ellipsoids  are  plotted  at  the  50%  probability  level.  The 
solvent  molecules  have  been  omitted  for  clarity. 
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known  E2  crystal  structures.  The  spread  of  over  12°,  with 
angles  found  distributed  over  the  entire  range,  indicate  the 
shallowness  of  the  potential  associated  with  the  bowing 
deformation.  In  a  system  as  dynamic  as  the  human  body,  the 
molecule  would  certainly  have  a  wide  range  of  conformations 
easily  available. 

The  packing  of  (I)  results  in  three-sided  cylinders  of 
estradiol  molecules  arranged  parallel  to  (101),  with  the  C18 
methyl  group  pointing  towards  the  center  (Fig.  2).  This,  of 
course,  aligns  the  large  hydrophobic  regions  of  the  molecules, 
as  well  as  the  hydrophilic  hydroxy  groups.  These  cylinders 
then  stack  on  top  of  each  other,  creating  three-sided  columns 
(Fig.  3)  held  together  by  hydrogen  bonds  involving  the 
hydroxy  groups  and  the  three  solvent  molecules.  A  more 
detailed  picture  of  the  hydrogen  bonding  between  the  solvent 
and  E2  molecules  can  be  seen  in  Fig.  4  and  from  the  data  in 
Table  3.  There  is  also  an  intercolumnar  hydrogen  bond,  which 
occurs  between  the  water  molecule  and  a  3-hydroxy  group. 
This,  along  with  the  intercolumnar  hydrogen-bonded  solvent, 


Figure  2 

The  molecular  arrangement  of  (I),  looking  down  the  three-sided  column. 
Displacement  ellipsoids  are  plotted  at  the  50%  probability  level 
[symmetry  codes:  (i)  x ,  y,  z;  (ii)  -2  -  x,  y  -  -1  -  z;  (iii)  x  -  1,  y,  z]. 


Figure  4 

The  hydrogen-bonded  region  of  the  column  shown  in  Fig.  3,  including  the 
solvent  molecules  and  only  the  A  or  D  rings  of  E2.  Displacement 
ellipsoids  are  plotted  at  the  50%  probability  level  [symmetry  codes:  (i)  x, 
y,  z;  (ii)  -x,y  - f,  -z;  (iii)  -l-x,y-\,  -z;  (iv)  1  +  x, y,  1  +  z;  (v)  -2  -  x, 
y  ~  -1  -  z;  (vi)  x,  y,  1  +  z;  (vii)  -1  -  X,  y,  z]. 
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Three-sided  column 


v/w*  Hydrophilic  layer 


Figure  5 

A  schematic  diagram  of  the  packed  columns  in  (I)  and  their  relation  to 
the  hydrophilic  layer  containing  the  hydroxy  groups  and  solvent 
molecules. 


creates  a  hydrophilic  layer  which  is  approximately  30°  from 
the  perpendicular  to  the  columns  (Fig.  5).  Structurally,  all 
bond  lengths  and  angles  are  in  expected  ranges  (Table  2). 


Figure  3 

The  packing  of  E2  groups  in  (I),  with  the  hydrogen-bonded  solvent 
molecules,  viewed  perpendicular  to  the  columns.  Dotted  lines  indicate  the 
region  detailed  in  Fig.  4.  Displacement  ellipsoids  are  plotted  at  the  50% 
probability  level  [symmetry  codes:  (i)  x ,  y,  z;  (ii)  1  +  x,  y,  1  +  z;  (iii)  -2  -  x, 
y  5*  1  z;  (iv)  -1  -  x,  y  -  *  -r,  (v)  x,  y,  1  +  z;  (vi)  -x,  y  -  \ ,  -z]. 


Experimental 

Crystals  were  grown  by  slow  evaporation  from  a  methanol  solution 
open  to  the  air.  It  is  presumed  that  the  methanol  was  either  wet  prior 
to  use  in  this  experiment  or  absorbed  atmospheric  moisture.  Crystals 
of  (I)  were  grown  over  the  course  of  3  d  and  harvested  when  reaching 
the  appropriate  size. 


Crystal  data 

Q8H2402  0.67CH40  0.33H20 
Mr  =  299.74 
Monoclinic,  P2X 
a  m  11.7152  (4)  A 
b  =  19.6270  (6)  A 
c  =  12.1310  (4)  A 
P  =  117.978  (1)° 

V  =  2463.34  (14)  A3 
Z  =  6 


Dx  =  1.212  Mg  nr3 
Ag  Ka  radiation 
Cell  parameters  from  7309 
reflections 
$  =  2.3-27.9° 

(i  =  0.05  mm"1 
T  =  298  (1)  K 
Cuboid,  colorless 
0.4  x  0.3  x  0.3  mm 
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Data  collection 

Siemens  SMART  Platform  CCD 
area-detector  diffractometer 
c o  scans 

Absorption  correction:  empirical 
via  multipole  expansion 
(Blessing,  1995)  using  SADABS 
(Sheldrick,  1996) 

T’min  =  0.840,  7max  =  0.985 

Refinement 

Refinement  on  F 2 
R[F2  >  ME2)]  =  0.058 
wR(F2)  m  0.145 
S  as  0.99 

21 122  reflections 
698  parameters 


44  195  measured  reflections 
21 122  independent  reflections 
15  597  reflections  with  I  >  2cr(7) 
Rim  -  0.048 
0max  «  27.9° 
h  =  —18  ->  19 
k  =  -32  29 

/  =  -20  ->  20 


H  atoms  treated  by  a  mixture  of 
independent  and  constrained 
refinement 

w  =  1  /[ffW)  +  (0.0812F)2] 
where  F  =  (F„2  +  2Fc2)/3 
(A/tr)m„  =  0.001 
APmax  =  0.41  e  A-3 
Apmi„  =  -0.23  e  A"3 


Table  1 

Comparison  of  ring-bowing  angles  (°)  in  reported  E2  crystal  structres. 


Compound 

Molecule  1 

Molecule  2 

Molecule  3 

angle 

angle 

angle 

17/8-estradiol-0.5H2Ot 

15.6 

17/J-estradiol-propanolt 

12.9 

17/3-estradiol-ureat 

5.6 

17y3-estradiol-0.5MeOHf 

10.4 

3.4 

0)5 

12.2 

11.9 

7.1 

t  Wiese  &  Brooks  (1994). 

$  Parrish  &  Pinkerton  (1999).  §  This  work. 

Table  2 

Selected  bond  lengths  (A). 

01A-C3A 

1.3716(18)  02 A - 

-cm 

1,4354  (19) 

01B—C3B 

1.3792  (17)  02 B- 

-072? 

1.4305  (19) 

01C-C3C 

1.3691  (18)  Q2C- 

-07  C 

1.4400  (19) 

Table  3 

Hydrogen-bonding  geometry  (A,  °). 

D-H-A 

D- H  H-  •  -A 

27..A 

D-H-.-A 

OL4— H1CXA-  •  *02/4 5 

0.74  (3)  1.89  (3) 

2.6226  (17) 

171  (3) 

OIB-HIOB-  .02^ 

0.85  (3)  1.83  (3) 

2.6654  (18) 

171  (3) 

OIC-HIOC -  03 

0.76(3)  1.82(3) 

2.5767  (18) 

170  (3) 

02A  — H2CX4  •  •  -042? 

0.76  (3)  1.97  (3) 

2.7017  (18) 

162  (2) 

02B~mOB  ■  -045" 

0.73  (2)  2.06  (2) 

2.7823  (19) 

168  (2) 

02C— H20C-  •  -OlC" 

0.73  (2)  1.99  (3) 

2.7169  (17) 

175  (3) 

03-H30A-  •  OMiv 

0.72  (4)  2.02  (4) 

2.743  (2) 

178  (4) 

03— H30B-*01£V 

0.76  (3)  2.04  (3) 

2.7926  (19) 

173  (3) 

04A-H40A  -  02C 

0.83(3)  1.85  (3) 

2.6809(18) 

174  (3) 

04B-H40B-  *  •04Avi 

0.78  (2)  1.88  (2) 

2.6462  (17) 

166  (2) 

Symmetry  codes:  (i)  1  +  x,  y,  1  +  z;  (ii)  -3  -  x,  \  +  y, 
~x,±  +  y,  ~z;  (v)  1  +  x,  y,  z\  (vi)  x  -  1 , y,  z. 


-2  —  z;  (iii)  x  -  1,  y,  z  —  1;  (iv) 


The  absolute  configuration  of  the  E2  molecules  in  (I)  was  known 
from  the  natural  product  starting  material.  The  intensity  data  were 
corrected  for  decay  and  absorption  using  SADABS  (Sheldrick,  1996). 
All  H  atoms  of  the  solvent  molecules  and  the  hydroxy  groups  of  the 
E2  molecules  were  located  in  a  difference  map  and  refined  with 
isotropic  displacement  parameters.  The  remaining  H  atoms  were 
included  with  idealized  geometries  (C-H  =  0.96-0.98  A),  and  their 
isotropic  displacement  parameters  were  refined. 

Data  collection:  SMART  (Bruker,  1998);  cell  refinement:  SAINT 
(Bruker,  1999);  data  reduction:  SAINT,  program(s)  used  to  solve 
structure:  SHELXS91  (Sheldrick,  1990);  program(s)  used  to  refine 
structure:  SHELXL91  (Sheldrick,  1997);  molecular  graphics: 
SHELXTL  (Siemens,  1994);  software  used  to  prepare  material  for 
publication:  SHELXTL . 

We  thank  the  College  of  Arts  and  Sciences  of  the  University 
of  Toledo  for  generous  financial  support  of  the  X-ray 
diffraction  facility. 


Supplementary  data  for  this  paper  are  available  from  the  IUCr  electronic 
archives  (Reference:  SQ1003).  Services  for  accessing  these  data  are 
described  at  the  back  of  the  journal. 
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Synopsis  The  initiation  of  a  comparative  charge  density  study  on  a  series  of  estrogen 
derivatives  revealed  the  need  for  a  standardized  local  coordinate  system  for  each  atom  in 
order  to  facilitate  comparison  of  the  multipole  refinements.  Herein  we  propose  such  a  system 
for  the  core  structure  of  a  generic  estrogen  molecule  and  suggest  starting  parameters  for  the 
population  of  the  most  prominent  multipoles. 

Abstract  A  comparative  charge  density  study  oh  a  series  of  estrogen  derivatives  has  been 
initiated.  The  study  utilizes  the  Hansen-Coppens  atom  centered  multipole  model  to  describe 
the  valence  electron  density  distribution.  Direct  comparison  of  the  population  parameters  for 
each  estrogen  after  the  respective  multipole  refinements  requires  standardization  of  the  atom 
centered  local  coordinate  systems.  Such  a  standard  coordinate  system  for  the  common 
estrogen  core  is  reported,  taking  advantage  of  the  shape  of  those  multipoles  which  have  the 
spatial  characteristics  of  sp2  and  sp3  hybrid  orbitals.  Additionally,  populating  these  principal 
multipoles  at  the  beginning  stage  of  the  refinements  improves  the  stability  of  these  large, 
highly  correlated  calculations. 

Keywords:  Estrogen;  charge  density;  coordinate  system 
1.  introduction 

Recent  advances  in  computers,  software,  and  X-ray  data  collection  technology  have 
made  charge  density  studies  on  a  series  of  larger  "small  molecules"  in  a  reasonable  amount  of 
time  feasible.  A  comparative  charge  density  analysis  has  been  initiated  to  study  a  series  of 
estrogen  derivatives. 

Electron  density  refinements  from  experimental  X-ray  diffraction  data  have  been 
performed  using  the  XD  suite  of  programs  (Koritsanszky  et  al.,  1995),  which  utilizes  the 
Hansen-Coppens  model  to  describe  the  valence  electron  density  (Coppens  1997,  Hansen  & 
Coppens  1978).  This  model  is  based  on  multipole  expansions  centred  on  the  atoms  and 


concomitantly  requires  a  user  defined  local  Cartesian  coordinate  system  for  each  atom  in  the 
structure . 

In  a  comparative  study,  consistency  in  this  model  setup  is  essential  to  enable  direct 
comparison  of  the  parameters  of  the  multipole  refinements  for  different  molecules  in  a  related 
series.  This  article  proposes  a  standard  local  coordinate  system  for  the  atoms  of  the  core 
structure  of  a  generic  estrogen  molecule. 

As  well  as  being  responsible  for  the  development  of  secondary  sexual  characteristics, 
it  has  been  determined  that  some  estrogens  are  responsible  for  the  initiation  and  progression 
of  certain  types  of  breast  cancer.  This  study  attempts  characterization  of  the  physical  and 
electronic  features  of  these  molecules  and  the  relationship  of  such  submolecular  properties  to 
biological  activity. 

Most  biologically  active  estrogens  are  A-ring  or  D-ring  isomers  of  the  naturally 
occuring  estradiol  or  derivatives  thereof,  containing  the  same  hydrocarbon  ring  core  structure. 
With  careful  choice  of  directions,  the  local  coordinate  system  can  take  advantage  of  the 
spatial  characteristics  of  certain  multipoles  that,  if  populated,  would  mimic  the  shape  of  the 
bonding  electron  density  derived  from  sp2  and  sp3  hybrid  orbitals..1  (Pichon-Pesme  et  al., 
1995;  Coppens,  1997)  This  allows  a  majority  of  the  valence  electron  density  to  be  modelled 
by  the  linear  combination  of  only  a  few  multipoles.  Given  the  number  of  parameters  for  such 
large,  non-centrosymmetric  systems,  typically  greater  than  40  atoms,  this  approach  can 
considerably  stabilize  the  initial  multipole  refinements.  By  using  a  standard  coordinate 
system,  upon  completion  of  the  multipole  refinement,  a  rational  comparison  can  then  be  made 
among  various  related  charge  density  studies. 

The  17^-estradiol  molecule,  along  with  the  labelling  scheme  of  the  core  structure,  is 
shown  in  scheme  1 .  For  all  methylene  groups,  hydrogen  atoms  on  the  same  side  of  the 
core  as  the  methyl  group  at  position  18  are  labelled  b,  those  on  the  other  side  are 
labelled  a.  The  C  and  D  rings  projected  down  the  C18-C13  bond  are  provided  as  an 
insert  to  define  the  orientation  of  the  methyl  group  hydrogen  atoms. 


1  It  is  important  to  distinguish  between  the  multipoles,  which  are  density  functions,  and  orbitals,  which 


are  derived  from  wavefunctions. 
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Scheme  1 


2.  Method 

The  cartesian  axes  x,y,z  of  the  above  mentioned  atom  centered  coordinate  systems  are 
calculated  in  the  program  XD  (Koritsanszky  et  al.,  1995)  based  on  two  user  defined  vectors. 
In  general,  the  first  vector  (vi)  for  a  given  atom  (1)  originates  at  its  nucleus,  is  directed 
towards  another  atom  (2)  in  the  system  -  or,  if  necessary,  a  "dummy"  atom  -  and  defines 
axis(l).  A  second  vector  (v2)  with  the  same  origin  points  towards  another  atom  (3),  or 
dummy  atom,  describing  the  axis(l)  -  axis(2)  plane.  A  third  vector  (v3)  is  taken  perpendicular 
to  this  plane, 

vi  =  (r2  -  n)  v2  =  (r3-ri)  v3=v,xv2 

where  1*1,  r2,  and  r3  are  the  position  vectors  of  atom(l),  atom(2),  and  atom(3),  respectively. 

An  orthonormal  vector  triplet  (eaxi,  e^,  e^)  is  then  formed  which  we  have  chosen  to  be  right 
handed. 

eaxi=Vi/|v1|  eax2=(v3XVi)/|(v3XVi)|  eaX3=v3/lv3| 

Although  axes  1, 2  and  3  can  be  assigned  labels  x,y,z  arbitrarily,  for  the  current  molecules  we 
define  axis(l)  as  x,  axis(2)  as  y,  and  axis(3)  as  z  for  non-hydrogen  atoms,  and  axis(l)  as  z, 
axis(2)  as  y  and  axis(3)  as  x  for  hydrogen  atoms. 

Given  these  assignments,  the  valence  electron  density  of  the  aromatic  sp2  type 
carbons  atoms  in  the  A-ring  (scheme  1),  can  be  initially  approximated  by  populating  two 
multipoles,  the  33+  and  20.  Pictorial  representations  of  the  angular  function  for  these 
multipoles  can  be  seen  in  Figure  1.  Defining  the  xy  plane  of  the  coordinate  system  to  lie  in 
the  plane  of  the  bound  atoms,  i.e.  the  aromatic  plane,  will  allow  the  positive  lobe  of  the  33+ 
octupole  to  align  with  the  three  covalent  bond  regions.  The  20  quadrupole  will  then  be 
positioned  perpendicular  to  this  plane,  a  negative  populations  adding  to  the  in  plane  o 
bonding  and  depopulating  the  p  orbital  as  required  for  ic  bond  formation. 


Figure  1  Coordinate  system  and  angular  functions  of  the  33+  and  20  multipoles  with  blue 
representing  positive  regions  and  red,  negative.  Functions  visualized  using  MATLAB  (2000). 

The  covalent  bonding  regions  of  the  sp3  type  carbon  atoms  in  the  B,  C,  and  D-rings  of 
the  generic  estrogen  can  effectively  be  modelled  by  a  single  multipole,  the  32-  octupole.  The 
axes  associated  with  the  32-  octupole  do  not  correspond  to  the  populated  regions,  hence  two 
"dummy"  atoms  must  be  used  to  define  the  local  coordinate  system  in  this  case.  Their 
positions  are  calculated  on  the  bisectors  of  the  bond  angles  (the  two-folds  axes  of  the 
tetrahedron)  at  the  atom  of  interest.  An  example  of  an  sp3  coordinate  system  and  the  32- 
octupole  is  shown  in  Figure  2. 

Approximation  of  the  bond  density  of  the  hydrogen  atoms  can  be  achieved  by 
initially  populating  a  single  dipole  along  z.  A  complete  list  of  vector  definitions  for  carbon 
and  hydrogen  atoms  for  the  estrogen  core  are  given  in  Tables  1  and  2.  Using  this  standardized 
set  of  coordinates  systems  for  future  theoretical  and  experimental  electron  density  studies  on 
estrogen  derivatives  (and  possibly  expanding  it  for  other  steroids  like  progesterone)  it  will  be 
possible  to  compare  individual  multipole  populations  for  different  molecules  in  this  series. 


Figure  2  Model  coordinate  system  for  a  sp3  type  atom  and  32-  multipole,  with  blue  representing 
positive  density  and  red,  negative.  (MATLAB,  2000). 

Populating  these  principal  multipoles  at  the  beginning  stage  of  the  refinements  may 
improve  the  stability  of  these  large,  highly  correlated  calculations.  We  have  thus  listed 
suitable  values  for  the  core  structure  based  on  our  own  refinements  (Chen  et  al.,  2003)  of  six 
different  estrogen  molecules  (Table  3).  We  note  that  the  sum  of  the  monopole  populations 
do  not  sum  to  zero  for  the  estrogen  core.  To  maintain  electroneutrality,  it  is  essential  that 
appropriate  starting  values  be  assigned  to  other  substituents,  e.g.  for  estradiol,  values  of -0.47 
and  0.29  would  be  appropriate  monopole  populations  for  the  O  and  H  atoms  of  the  hydroxy 
groups. 

3.  Conclusion 

Examples  of  typical  input  for  the  XD  program  and  pictorial  representation  of  the 
coordinate  schemes  for  the  hydrocarbon  core  of  estrogens  as  recommended  in  this  paper  are 
available  as  supplementary  material. 


Table  1  Atoms  used  to  define  vectors  Vi  and  v2  for  carbon  atoms3 


Atom(l) 

Atom  (2) 

Atom(3) 

Cl 

C2 

CIO 

C2 

C3 

Cl 

C3 

C4 

C2 

C4 

C5 

C3 

C5 

CIO 

C4 

C6 

Duml  (C5  -  C7) 

Dum2  (C5  -  H6a) 

Cl 

Dum3  (C6  -  C8) 

Dum4  (C6  -  H7a) 

C8 

Dum5  (Cl  -  C9) 

Dum6  (C7  -  C14) 

C9 

Dum7  (C8  -  Cll) 

Dum8  (C8  -  CIO) 

CIO 

Cl 

C5 

Cll 

Dum9  (C9  -  C12) 

DumlO(C12-Hlla) 

C12 

Dumll  (Cll -C13) 

Duml 2  (C13  -  H12a) 

C13 

Duml3  (C12-  C14) 

Duml4  (C14-C17) 

C14 

Duml5(C8-C15) 

Duml6  (C8  -  C13) 

C15 

Duml7  (C14  -  C16) 

Duml8  (C14  -  H15a) 

C16 

Duml  9  (C15  -  C17) 

Dum20  (C15  -  H16a) 

C17 

Dum21  (C13-C16) 

Dum22  (C16  -  02) 

C18 

Dum23  (H18a  -  H18b) 

Dum24  (HI  8b  -  H18c) 

a  Bond  lengths  of  atoms  used  to  calculate  dummy  atom  positions  (in  parentheses)  are  adjusted  to  be 
equidistant  from  the  atom  of  interest  to  ensure  a  true  bisection. 


Table  2.  Atoms  used  to  define  vectors  Vj  and  v2  for  hydrogen  atoms 


Atom(l) 

Atom(2) 

Atom(3) 

HI 

Cl 

C9 

H2 

C2 

HI 

H4 

C4 

01 

H6a 

C6 

H6b 

H6b 

C6 

H6a 

H7a 

Cl 

H7b 

H7b 

Cl 

H7a 

H8 

C8 

03 

H9 

C9 

H14 

Hlla 

Cll 

Hllb 

Hllb 

Cll 

Hlla 

H12a 

C12 

H12b 

H12b 

C12 

H12a 

H14 

C14 

H9 

H15a 

C15 

H15b 

H15b 

C15 

HI  5a 

H16a 

C16 

H16b 

H16b 

C16 

H16a 

HI  7 

C17 

02 

H18a 

C18 

H18b 

H18b 

C18 

H18c 

H18c 

C18 

HI  8a 

Table  3:  Initial  set  of  population  parameters  for  the  most  important  multipoles. 


spi 

sp3 

Atom 

00 

33+ 

32- 

00 

10 

Cl 

-0.22(3) 

0.32(5) 

0.24(4) 

0.14(4) 

C2 

-0.28(7) 

-0.21(5) 

0.34(5) 

H2 

0.24(4) 

0.15(4) 

C3 

0.34(7) 

mm 

C4 

0.33(4) 

u 

-0.13(3) 

H7x 

0.18(4) 

0.14(2) 

mm 

-0.27(4) 

H8 

0.19(4) 

0.13(5) 

C7 

-0.28(6) 

0.31(6) 

H9 

0.18(2) 

0.12(4) 

-0.18(6) 

0.38(5) 

Hllx 

0.18(4) 

-0.16(6) 

0.33(9) 

H12x 

0.18(4) 

CIO 

0.36(4) 

Cll 

C12 

-0.28(5) 

0.32(6) 

H16x 

0.17(5) 

0.15(3) 

C13 

-0.21(6) 

H17 

0.10(2) 

0.18(2) 

C14 

-0.15(4) 

0.14(3) 

0.13(4) 

C15 

-0.32(7) 

0.30(7) 

C16 

-0.35(5) 

0.32(7) 

IHIH 

0.16(4) 

0.30(5) 

C18 

-0.39(6) 

0.22(6) 

Hydrogen  labels  such  as  H6x,  represent  both  H6A  and  H6B.  Populations  are  average  values  from  6 
studies  on  *  np-estradiol.V^MeOH,  17p-estradiol.urea,  17a-estradiol,  17a-estradiol.1/2H20,  16a,  17p  - 
estriol,  estrone,  k’s  for  H  atoms  were  fixed  at  1.4.  Other  values  for  k  would  be  expected  to  introduce 
small  changes  to  the  refined  monopole  populations,  but  to  have  minimal  effect  on  the  higher  poles. 
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In  order  to  relate  molecular  electrostatic  potential  to  biological  activities  of  estrogens,  a 
comparative  charge  density  study  of  different  derivatives  has  been  initiated.  The  first  completed 
charge  density  analysis  of  this  series  for  17.y9-estradiol«urea  is  presented  here.  It  has  been 
demonstrated  that  it  is  possible  to  perform  a  quality  experimental  charge  density  study  for  the  big 
organic  system  with  56  atoms  in  the  non-centrosymmetric  space  group.  The  special  tools  such  as 
the  optimal  coordinate  system  and  slow  initially  constrained  refinement  have  been  used  to 
accomplish  this  study. 

Our  results  for  the  urea  molecule  reasonably  agree  with  the  previous  experimental  and 
theoretical  results  showing  the  quality  of  our  data.  In  the  17/?-estradiol  molecule,  oxygen  atoms 
appeared  to  be  sp3  in  shape,  exhibiting  two  consistent,  distinct  lone  pairs  despite  of  different 
chemical  environments.  No  interaction  of  hydroxy  group  oxygen  lone  pairs  with  the  #  orbitals  of 
the  aromatic  ring  has  been  observed.  The  analysis  of  the  electrostatic  potential  revealed  that  the 
negative  potential  in  the  lone  pair  region  of  the  oxygen  atoms  is  quite  different.  The  Bader's 
topological  analysis  has  been  performed,  and  atomic  charges  have  been  estimated.  The  results 
are  compared  and  discussed.  ' 

Estrogens  are  known  to  be  responsible  for  the  development  of  secondary  sexual  characteristics, 
as  well  as  they  effect  growth,  differentiation  and  function  of  wide  range  of  tissues  [1].  It  has  also 
been  shown  [2]  that  some  estrogens  are  responsible  for  the  initiation  and  progression  of  certain 
types  of  breast  cancer.  These  molecules  have  the  ability  to  bind  as  ligands  to  the  estrogen 
receptor  in  the  first  of  many  steps  which  could  result  in  the  activation  (agonistic  effect)  or 
repression  (antagonistic  effect)  of  genes  critical  in  the  mechanism  of  tumor  growth.  While  the 
mechanism  by  which  the  estrogens  influence  cancer  is  currently  unknown,  subtle  changes  in  the 
chemical  structure  of  estradiol  and  the  other  estrogens  are  known  to  elicit  different  biological 
responses  in  the  development  of  cancer.  It  has  been  suggested  [3]  that  the  agonistic/antagonistic 
responses  of  the  different  estrogens  can  be  related  to  such  physical  properties  as  their 
electrostatic  potential  and  charge  density  distribution. 

A  comparative  charge  density  study  of  a  series  of  estrogen  derivatives  has  been  initiated  in  our 
group,  with  the  primary  goal  of  relating  the  electronic  and  physical  structure  of  the  molecule  to 
biological  action.  In  addition  to  molecular  properties,  charge  density  studies  provide  significant 
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fundamental  information  about  the  functional  groups  and  the  atoms  which  make  them  up, 
helping  to  investigate  basic  chemical  principles  of  small  organic  molecules  in  general.  This  paper 
will  pay  particular  attention  to  the  hydroxy  groups  of  the  estrogen  molecule.  The  hydroxy  groups 
are  some  of  the  most  chemically  interesting  portion  of  the  estrogen  molecule,  and  their 
interaction  with  the  receptor  is  thought  to  be  mostly  responsible  for  their  activity. 

Described  here  is  the  charge  density  study  of  17/?-estradiol»urea.  This  particular  crystal  system 
was  chosen  to  be  the  first  in  the  series  for  two  reasons.  First,  17/?-estradiol  (E2)  is  the  most 
common  naturally  occurring  estrogen  and  is  typically  the  standard  by  which  all  other  estrogenic 
activities  are  related.  E2  demonstrates  an  agonistic  response  in  in  vitro  tests  which  typically 
utilize  the  estrogen  receptor  in  either  competitive  yeast  or  human  cell  assays  [4-5].  Secondly,  it 
contained  a  solvated  urea  molecule  which  has  been  extensively  studied.  Both  theoretical  and 
experimental  charge  density  studies  have  been  performed  on  this  molecule  [6-7].  This  provides  a 
direct  means  of  comparison  to  test  our  methods  against  previous  results. 

Experimental 

Data  Collection  and  Reduction 

Crystals  were  grown  by  slow  evaporation  as  described  by  Duax  in  the  original  publication  of 
the  crystal  structure  [8].  The  clear,  colorless  crystal  of  approximate  dimensions  0.35  x  0.37  x 
0.40  mm  was  attached  to  a  30  //m  carbon  fiber  using  a  small  amount  of  epoxy  resin.  Data  were 
collected  at  100.0(1)  K  on  a  Bruker  Platform  Diffractometer  equipped  with  a  SMART  6000 
CCD  area  detector  located  at  7.12  cm  from  the  goniometer  center.  The  Oxford  Cryostream 
system  has  been  used  to  cool  the  crystal.  The  crystals  were  irradiated  with  graphite 
monochromated  Mo  Ka  radiation  (X  =  0.7107  A).  0.3°  omega  scans  were  carried  out  at  four  phi 
settings  and  three  detector  positions  in  20,  resulting  in  more  than  7400  frames.  A  frame  time  of 
60  seconds  was  used  for  the  low  angle  setting,  and  180  seconds  for  the  medium  and  high  angle 
settings.  A  complete  experimental  protocol  is  deposited. 

The  intensities  were  integrated  using  the  program  SAINT  [9].  After  several  integrations  and 
preliminary  refinements,  it  was  found  that  different  integration  parameters  were  required  for 
each  detector  setting  to  obtain  the  best  possible  structure  factors.  The  box  sizes,  profile  fitting 
parameters,  and  the  simple  sum  perimeter  limits  were  adjusted  for  each  detector  setting.  The  unit 
cell  parameters  have  been  finally  refined  in  SAINT  using  a  total  of  9999  reflections. 

Lorentz-polarization  and  A/2  corrections  were  applied  to  the  integrated  intensities  with  a  A/2 
correction  factor  of  0.0012,  this  value  being  taken  from  the  earlier  study  [10].  The  program 
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SORTAV  [11]  was  used  for  outlier  determination,  averaging,  and  applying  the  absorption 
correction.  Only  reflections  with  I  >  3o(I)  and  below  1.00  A'1  were  used  during  the  multipole 
refinements. 


Table  1.  The  crystal  information  and  experimental  details. 


Formula 

C19H28N2O3 

Formula  Weight 

332.43 

Crystal  System 

Orthorhombic 

Space  Group 

P2,2i2i 

a 

7.9022(9)  A 

b 

9.2228(10)  A 

C 

24.5890(28)  A 

a=p=y 

90.0° 

Volume 

1792.06  A3 

Z 

4 

T 

100.0  K 

X 

0.7107  A 

(sin  $^)max 

1.180  A’1 

Reflections  Collected 

110999 

Rejected  Outliers 

779 

Unique  Reflections 

13187 

Included  in  the  Refinement  [/>  3  <7(7)] 

7626 

Data  completeness 

98.6  % 

Average  redundancy 

8.4 

Rint 

0.0411 

Ri  (spherical) 

0.065 

WR2  (spherical) 

0.145 

R(F2)  [/>  3 o(I)]  (multipole) 

0.029 

R(F2)  all  data  (multipole) 

0.036 

wR  (F2) 

0.045 

Weighting  scheme 

I/O2 

Goodness  of  Fit 

1.393 

Least-Squares  Refinements 

The  crystal  structure  was  re-solved,  and  the  spherical  atom  model  was  refined  using  the 
program  SHELXL-97  [12].  17P-estradiol  and  urea  molecules  are  shown  in  Figure  1.  All 
hydrogen  atoms  were  located  from  the  difference  Fourier  map.  The  results  of  this  refinement 
served  as  the  starting  point  for  the  aspherical  atom  refinements  that  were  performed  using  the  XD 
software  package  [13]. 

C118I 


0(3) 


Figure  1.  17P-estradiol  and  urea  molecules  showing  50%  probability  thermal  ellipsoids  as  at 
100K.  Hydrogen  atoms  are  omitted  for  clarity. 

The  aspherical  atom  model  [14],  used  in  further  refinements  consists  of  three  components.  The 
first  and  second  describe  spherically  averaged  core  and  valence  densities  normalized  to  one 
electron,  while  the  third  one  describes  the  aspherical  valence  density.  p(r)  is  the  electron  density, 
Pa  Pv,  and  Pim  are  population  coefficients,  the  core  population  (Pc)  was  fixed  in  this  study. 
Asphericity  is  taken  into  account  by  expanding  of  the  electron  density  over  real  combinations  of 
spherical  harmonics,  yim,  which  satisfy  local  symmetry.  The  variable  radial  parameters,  rand  x', 
allow  for  expansion  or  contraction  of  the  valence  terms.  Ri  are  normalized  Slater-type  radial 
functions. 

Patomic(r)  —  PcPcore  +  PvK  3  p  valence  (Kr) 

4  1 

+  ^  k'  3  Ri(k'  r )  ^  Pimyim(r  /  r ) 

1=1  m=-l 
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Because  the  model  is  based  on  atom  centered  multipoles,  a  coordinate  system  is  required 
for  each  atom  in  the  system.  The  core  structure  of  the  estrogen  molecule  was  set  up  using  the 
standard  estrogen  coordinate  system  to  take  advantage  of  the  shape  of  certain  multipoles  to 
increase  the  stability  of  the  initial  multipole  refinements  [15].  The  vectors  for  the  oxygen  atoms 
of  the  estrogen  molecule  were  directed  to  bound  atoms.  The  coordinate  system  setup  for  the  urea 
molecule  is  deposited.  - 

Before  the  refinements  were  begun,  the  hydrogen  atom  positions  were  corrected  by  extending 
the  X  -  H  distances  to  their  neutron  bond  lengths  (Csp2H  =  1.08  A,  Csp3H  =  1.10  A,  Csp3H2  =  1.09 
A,  Csp3H3  =1.06  A,  OH  =  0.97  A,  Nsp2H2  =1.01  A).  These  values  were  fixed  in  the  refinement. 
A  well  known  correlation  exists  between  the  thermal  parameter  and  the  kappa  parameters  of  a 
hydrogen  atom,  only  allowing  the  kappa  parameters  to  be  refined  when  the  thermal  parameters 
of  the  hydrogen  atom  are  known  [16].  To  avoid  this  problem,  fixed  kappa  values  (kh  =  1.20  and 
X#  =1.20)  were  used,  while  the  thermal  parameters  were  allowed  to  refine  [17-18]. 

Initially,  to  reduce  the  number  of  parameters  in  the  least  squares  process,  chemical  constraints 
and  only  the  most  significant  multipoles  were  considered.  By  the  end  of  the  refinement,  all  non¬ 
hydrogen  atoms  were  treated  up  to  the  octupole  level  with  no  chemical  constraints.  The 
hydrogen  atoms  were  refined  with  only  two  multipoles  having  cylindrical  symmetry  along  the  X 
-  H  bond.  Hydrogen  atoms  bound  to  the  same  carbon  have  been  constrained  to  be  equal.  The 
electroneutrality  constraint  has  been  imposed  on  both  of  two  molecules.  Complete  tables  giving 
the  final  structural  information,  such  as  atomic  coordinates,  Uy’s,  IWs,  bond  lengths,  bond 
angles  and  multipole  parameters  are  deposited. 

Topological  Analysis,  Atomic  Charges  and  Electrostatic  Potential 

The  topological  analysis  was  performed  based  on  the  Atoms  in  Molecules  (AIM)  Quantum 
Theory  [19].  The  XD  property  assessment  program,  XDPROP,  was  used  to  perform  a 
quantitative  critical  point*  analysis  on  all  covalent  and  hydrogen  bonding  interactions,  as  well  as 
oxygen  lone  pair  densities.  The  atomic  basin  boundaries  were  defined  and  the  atomic  charges 
(Bader's  charges)  were  integrated  within  the  basins  with  the  program  TOPXD  [20],  The 
integrated  atomic  charge  q(£2)  over  the  atomic  basin  Q,  is  defined  as  a  difference  between  nuclear 
Zn  and  electronic  N(Gl)  charges: 


*  Critical  points  are  special  points  in  the  electron  density,  where  Vp(rcp)=0.  The  saddle  (3,-1)  point  reflects  the 
chemical  bonding.  In  the  Laplacian,  the  (3, +3)  minimum  point  corresponds  to  the  maximum  curvature  of  the 
electron  density  in  the  charge  concentration  regions.  See  ref.  [19]. 
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q(Q)  =  Zq  -  N(Q),  N( Q)=J  p(r)dx. 

a 

This  method  allows  one  an  exhaustive  partitioning  of  the  molecules  into  the  respective  atoms, 
leading  to  the  additivity  and  transferability  of  it's  properties  [21]. 

The  calculation  of  the  electrostatic  potential  (ESP)  was  performed  using  the  program 
XDPROP.  It  can  be  calculated  by  the  evaluation  of  \  . 


V(r)  =  £ 

A 


ZA 

RA-r| 


fP(r')dr' 

J  |r'_r 


where  ZA  is  the  charge  on  nucleus  A  located  at  RA.  The  construction  assembles  positive  charges 
at  the  nuclei  and  a  continuous  distribution  of  electrons.  The  calculations  were  performed  based 
on  the  isolated  from  the  crystal  P-estradiol  molecule. 


Results  and  Discussion 

Packing 

The  molecular  structure  and  hydrogen  bonding  pattern  were  found  to  be  the  same  as  described 
by  Duax  [8],  The  packing  of  this  structure  results  in  a  formation  of  six  hydrogen  bonds  ranging 
in  donor-acceptor  length  from  2.64  -  3.10  A,  these  are  listed  in  Table  2.  The  0(1)  atom  acts  as 
an  acceptor  in  two  hydrogen  bonds  with  hydrogens  of  two  independent  urea  molecules  (H(1NB) 
and  H(2NB)).  The  HO  vectors  of  these  hydrogen  bonds  are  pointed  almost  directly  at  the 
positions  of  the  lone  pairs  of  0(1).  The  0(2)  atom  has  a  very  similar  hydrogen  bonding  structure 
to  0(1)  in  which  it  acts  as  an  acceptor  in  two  hydrogen  bonds  and  a  donator  in  one.  The 
hydrogen  atoms  from  the  0(1)  atom  of  another  independent  E2  molecule  and  H(2NB)  of  a  urea 
molecule  line  up  almost  directly  with  the  lone  pair  positions  of  0(2).  Topological  analysis 
yielded  in  a  (3,-1)  critical  point  and  a  bond  path1  for  each  hydrogen  -  acceptor  interaction, 
verifying  the  true  bonding  interaction  in  these  hydrogen  bonds  [22].  The  complete  table  of 
results  from  hydrogen  bond  topology  is  deposited.  . 


'  Bond  path  is  a  line  along  which  the  electron  density  decreases  for  any  lateral  displacement. 
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Table  2.  The  geometrical  information  for  the  hydrogen  bonds. 


D— H-A 

d(D  -A),  A 

Z  DHA,  ° 

0(2)— H(20)  •  *0(3) 

2.6395(7) 

169.1(16) 

0(1)-H(10)-  -0(2) 

2.6698(7) 

171.0(20) 

N(l)— H(1NA)-  -0(3) 

2.8231(10) 

168.6(18) 

N(l)— H(1NB)-  *0(1) 

2.9783(9) 

135.2(16) 

N(2)— H(2NB)  •  •  0(1) 

3.0086(9) 

161.1(17) 

N(2)— H(2NA)  •  •  0(2) 

3.1023(8) 

169.1(16) 

Kappa  Parameters 

The  coefficients  of  expansion/contraction  of  the  spherical  valence  and  deformation 
densities,  K  and  k’  parameters  respectively,  are  shown  in  Table  3.  The  grouping  of  individual 
atoms  were  based  on  the  atom  type,  hybridization  and  chemical  environment.  A  total  of  twelve 
kappa  sets  were  utilized  in  order  to  allow  the  necessary  flexibility  while  attempting  to  maintain  a 
minimum  number  of  parameters.  As  mentioned  earlier,  the  hydrogen  atom  kappa  values  were 
fixed  at  1.20.  All  other  kappa  values  were  allowed  to  refine  during  the  least  squares 
minimization. 

Table  3.  The  expansion/contraction  (k  and  K')parameters  for  the  17p-estradiol«urea  crystal 


Atoms 

Kappa 

K 

K' 

0(1),  0(2) 

1 

0.971 

0.854 

C(3) 

2 

1.009 

0.879 

C(17) 

3 

1.012 

0.933 

C(l),  C(2),  C(4) 

4 

0.976 

0.896 

C(5),  C(10) 

5 

0.986 

0.838 

C(6),  C(7),  C(8),  C(9), 

C(ll),  C(12),  C(13),  C(14), 

6 

0.985 

0.878 

C(15),  C(16),  C(17),  C(18) 

all  C  -  H  hydrogen  atoms 

7 

1,200 

1.200 

8 


H(10),  H(20) 

8 

1.200 

1.200 

C(19) 

9 

0.986 

0.818 

0(3) 

10 

0.991 

1.080 

N(1),N(2) 

11 

0.996 

1.004 

all  N  -  H  hydrogen  atoms 

12 

1.200 

1.200 

Atomic  Charges 

The  net  atomic  charges  were  derived  directly  from  the  monopole  populations  of  the  multipole 
model,  as  well  as  from  integrating  of  the  electron  density  within  the  atomic  basins  as  described 
by  the  AIM  theory.  A  complete  list  of  atomic  charges  determined  by  both  methods  are  listed  in 
Table  4. 

As  mentioned  earlier,  the  urea  molecule  gives  us  the  opportunity  to  compare  our  results  with 
previous  theoretical  and  experimental  results  [6-7].  The  atomic  charges  derived  directly  from  the 
monopole  populations  of  the  multipole  model  compared  favorably  to  the  experimental  values 
reported*  by  Zavodnik  et  al.  [6],  with  values  matching  within  one  e.s.d.  Bader’s  charges  are  in  a 
fair  agreement  with  those  reported  by  Gatti  et  al.  [7].  Generally,  Bader's  charges  are  significantly 
bigger  than  the  charges  derived  from  the  monopole  populations;  this  is  not  surprising  due  to  the 
different  ways  of  calculations.  Upon  first  inspection,  a  charge  of  +2.13  e'  for  the  central  carbon 
seems  tremendously  high.  At  the  same  time,  there  is  evidence  in  the  literature,  from  both 
theoretical  and  experimental  data,  that  this  is  not  unreasonable  for  a  carbon  surrounded  by 
electronegative  oxygen  and  nitrogen  atoms  [7,  23]. 

The  net  atomic  charges  for  the  atoms  of  the  estrogen  molecule  derived  from  the  monopole 
populations  of  the  multipole  model  were  as  expected  for  the  given  atom  types  and  their  chemical 
neighbors.  In  addition,  like  atoms  with  similar  chemical  environments  exhibited  very  consistent 
monopole  populations.  Again,  the  charge  integrations  over  atomic  basins  yielded  much  greater 
charges  for  the  oxygen  atoms  and  the  atoms  directly  bound  to  them,  maintaining  the  consistency 
between  chemically  equivalent  atoms. 

Table  4.  Net  atomic  charges  for  the  atoms  of  17(3-estradiol  •  urea 


*  See  the  deposition  material  of  Zavodnik  et  al.  (1999)  Acta  Crystallogr.  B55,  45-54. 
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Atom 

q(Pv) 

q(Q) 

Atom 

q(Pv) 

q(O) 

0(1) 

-0.534(9) 

-1.24 

H(10) 

0.401(9) 

0.68 

0(2) 

-0.527(9) 

-1.11 

H(20) 

0.391(9) 

0.67 

C(l) 

-0.233(16) 

-0.19 

H(l) 

0.206(10) 

0.14 

0(2) 

-0.227(16) 

-0.12 

H(2) 

0.248(9) 

0.20 

C(3) 

0.148(14) 

0.49 

H(4) 

0.249(9) 

0.21 

0(4) 

-0.267(15) 

-0.21 

H(6A) 

0.147(7) 

0.15 

0(5) 

-0.110(15) 

-0.09 

H(6B) 

0.147(7) 

0.16 

0(6) 

-0.265(16) 

-0.28 

H(7A) 

0.146(7) 

0.08 

0(7) 

-0.279(16) 

-0.15 

H(7B) 

0.146(7) 

0.10 

0(8) 

-0.121(14) 

-0.12 

H(8) 

0.180(9) 

0.16 

0(9) 

-0.131(15) 

-0.07 

H(9) 

0.180(9) 

0.09 

0(10) 

-0.106(15) 

-0.05 

H(11A) 

0.147(7) 

0.08 

0(11) 

-0.277(15) 

-0.16 

H(11B) 

0.147(7) 

0.09 

0(12) 

-0.273(15) 

-0.14 

H(12A) 

0.141(7) 

0.07 

0(13) 

-0.158(16) 

-0.08 

H(12B) 

0.141(7) 

0.06 

0(14) 

-0.111(14) 

-0.04 

H(14) 

0.176(8) 

0.11 

0(15) 

-0.294(15) 

-0.23 

H(15A) 

0.141(7) 

0.10 

0(16) 

-0.273(16) 

-0.18 

H(15B) 

0.141(7) 

0.13 

0(17) 

0.168(13) 

0.36 

H(16A) 

0.146(7) 

0.11 

0(18) 

-0.402(17) 

-0.24 

H(16B) 

0.146(7) 

0.14 

0(3) 

-0.234(9) 

-1.27 

H(17) 

0.120(9) 

0.14 

N(l) 

-0.282(12) 

-1.41 

H(18A) 

0.122(6) 

0.06 

N(2) 

-0.286(11) 

-1.40 

H(18B) 

0.122(6) 

0.06 

0(19) 

-0.030(14) 

2.13 

H(18C) 

0.122(6) 

0.06 

H(1NA) 

0.205(10) 

0.49 

H(1NB) 

0.209(11) 

0.49 

H(2NA) 

0.207(10) 

0.47 

H(2NB) 

0.208(10) 

0.49 

q(Pv)  -  net  atomic  charges  derived  directly  from  the  monopole  populations  of  the  multipole  model; 
q(£2)  -  Bader's  net  atomic  charges 


Topological  Analysis 
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a)  urea  molecule 

Overall  values  of  the  electron  density  at  the  bond  critical  points  are  in  a  fair  agreement  with 
those  reported  by  Zavodnik  et  al.  [6]  and  Gatti  et  al.  [7],  with  the  exception  of  the  N  -  H  bonds 
(Table  5).  Zavodnik  et  al.  [6]  reported  significant  differences  (~0.5  eA'3)  between  the  theoretical 
and  experimental  electron  densities  at  the  (3,-1)  critical  point  for  the  N  -  H  bond  of  the  urea 
molecule.  The  pi r)  values  in  this  study  are  much  closer  to  the  theoretical  values.  The  difference 
in  the  Laplacian  values  at  (3,-1)  critical  points  between  experimental  and  theoretical  values  is 
much  bigger,  showing  the  sensitivity  of  the  electron  density  curvature  to  the  different  approaches 
of  it's  estimation,  including  the  deficiency  of  the  multipole  model  [24]. 


Table  5.  Topological  analysis  for  the  urea  molecule. 


this  work  Zavodnik,  et  al.  [6]  Gatti,  et  al.  [7] 

experiment  experiment  theory  theory 


bond 

C-0 


(C-N)average 

(N-H)average 


pir) 

pi r) 

Pit) 

2.79 

-33.15 

2.54 

-18.86 

2.57 

-7.92 

2.57 

-7.95 

2.26 

-24.20 

2.54 

-37.66 

2.35 

-27.61 

2.36 

-27.71 

2.22 

-28.91 

1.82 

-31.91 

2.34 

-47.06 

2.34 

-47.23 

The  units  for  the  electron  density,  p{ r),  are  eA'3,  and  for  the  Laplacian,  V2p(r),  they  are  eA'5. 


A  critical  point  search  was  also  performed  in  the  Laplacian  to  find  the  charge  density 
concentration  associated  with  the  lone  pairs  on  the  urea  oxygen.  The  critical  points  are  (3, +3)  in 
type  indicating  that  V2pir)  is  a  local  minimum  at  rc.  The  critical  point  information,  as  well  as  a 
plot  of  the  Laplacian  in  the  plane  of  the  lone  pairs  of  the  6(3)  atom  can  be  seen  Table  6  and  Fig. 
2  respectively.  The  angle  between  the  lone  pairs,  being  greater  that  140°,  is  much  larger  than  the 
ideal  trigonal  planar  geometry  in  agreement  with  VESPR  theory.  The  inconsistency  between  the 
C(3)  -  0(1)  -  LP  angles  are  most  likely  due  to  the  shallow  curved  shape  of  the  Laplacian  in  the 
lone  pair  region.  An  extremely  small  change  in  V2p( r)  could  easily  move  the  critical  point  by  5  - 
10°.  Gatti  et  al.  [7]  also  found  two  (3, +3)  critical  points  in  the  Laplacian,  corresponding  to  the 
urea  oxygen  lone  pairs.  Although  one  of  their  points  is  very  similar  to  those  in  Table  6,  having 
the  V2pc( r)  =  -137  eA'5  and  the  0(1)-LP  distance  of  -0.339  A,  another  (3, +3)  point  is 
significantly  away  from  the  oxygen  nucleus  (-0.383  A)  having  much  lower  value  of  the 
Laplacian  (~  -48  eA'5).  As  it  has  been  shown  before,  the  Laplacian  values  (and  it's  curvatures) 
are  very  'method-sensitive'. 
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Table  6.  Lone  pair  Laplacian  critical  point  information  for  the  0(3)  atom. 


Atom 

V2p(rc) 

(eA'5) 

CP 

Type 

0(1)-LP 

(A) 

LP1-0(1)-LP2 

(°) 

C3-0(1)-LP 

(°) 

0(3) 

LP1 

-149 

(3, +3) 

144 

118 

LP2 

-148 

(3, +3) 

144 

92 

Figure  2.  The  negative  part  of  the  Laplacian  in  the  plane  of  the  0(3)  atom  lone  pairs.  Contour 
interval  is  -5  eA'5  starting  at  -100  eA'5. 

b)  E2  molecule 

The  critical  point  analysis  of  covalent  bonds  of  the  estrogen  molecule  demonstrated  the 
characteristic  (3,-1)  critical  points,  with  expected  p(rc),  V2p(rc),  and  ellipticity  values  given  the 
hybridization  and  connectivity.  The  complete  list  of  all  bond  critical  point  properties  is 
deposited.  The  average  p( rc)  values  are  1.91  eA'3  for  the  C  -  O  bond,  2.17  eA'3  for  Car  -  Car,  1-66 
eA 3  for  C  —  C,  1.91  eA 3  for  C  —  H  and  2.07  eA 3  for  the  O  —  H  bond.  The  critical  points  are  well 
centered  in  the  homonuclear  C  -  C  bonds,  and  proportionally  displaced  away  from  the  more 
electronegative  atom  in  the  heteronuclear  bonds. 

Critical  points  in  the  Laplacian  corresponding  to  the  lone  pairs  of  both  hydroxy  oxygen  atoms 
are  of  the  special  interest.  Table  7  indicates  the  positions  of  the  lone  pair  critical  points  in  the 
Laplacian.  The  deformation  electron  density  map  of  the  lone  pairs  of  the  oxygen  0(1),  bound  to 
the  aromatic  ring  is  presented  in  Fig.  3,  while  plot  of  the  Laplacian  in  the  plane  of  the  lone  pairs 
of  the  0(2)  atom  is  shown  in  Fig.  4.  The  finding  of  two  distinct  lone  pair  regions  after  the 
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multipole  refinement  for  0(1)  was  somewhat  surprising.  We  expected  that  the  aromatic  hydroxy 
group,  in  which  hydrogen  lies  in  the  plane  of  the  aromatic  ring,  would  be  of  sp2  type  geometry 
having  some  interaction  with  the  JV-  orbitals  of  the  aromatic  ring  (the  same  as  in  phenol), 
however,  the  results  of  this  study  indicate  that  this  is  not  the  ground  state  configuration  in  the 
estrogen  molecule.  The  distances  of  the  lone  pairs  from  the  oxygen  atoms  are  very  consistent. 
This  is  believed  to  be  at  least  partially  to  due  to  the  fact  that  the  contraction/expansion 
coefficients  for  both  oxygen  atoms  are  constrained  to  be  equal.  Same  as  in  the  urea  molecule, 
the  angles  between  the  lone  pairs  agree  with  VESPR  theory  in  that  they  are  much  larger  than  the 
ideal  tetrahedral  value. 


Table  7.  Lone  pair  Laplacian  critical  point  information  for  the  0(1)  and  0(2)  atoms. 


Figure  4.  The  negative  part  of  the  Laplacian  in  the  plane  of  the  lone  pairs  of  the  0(2)  atom. 
Contour  intervals  are  -5  eA'5  starting  at  -80  eA'5. 


Electrostatic  Potential 

The  core  structure  of  the  estrogen  molecule  generally  is  hydrophobic  and  relatively  predictable 
in  terms  of  die  electrostatic  potential.  The  ESP  plotted  on  the  van  der  Waals  molecular  surface 
for  the  E2  molecule  with  the  program  package  MOLMOL  [25]  is  shown  in  Fig.  6.  Since  the 
isolated  molecule  (taken  from  the  crystal)  is  calculated  in  this  case,  the  ESP  does  not  reflect  the 
interactions  in  the  condensed  phase,  such  as  hydrogen  bonding.  The  negative  potential  areas 
(red)  are  produced  by  the  aromatic  ring,  oxygens  and  some  of  carbon  atoms  (Fig.  6),  while  the 
hydrogen  atoms  involved  in  the  hydrogen  bonding  produce  highly  positive  areas  (blue)  on  the 
molecular  surface. 


Figure  6.  The  E2  electrostatic  potential  plotted  on  the  van  der  Waals  molecular  surface.  Negative 
values  are  in  red,  positive  are  in  blue. 


The  hydroxy  groups  are  of  particular  interest  in  the  E2  molecule,  offering  varying 
environments  being  bound  to  sp2  and  sp3  type  carbons  and  being  involved  in  several  hydrogen 
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bonds.  The  consistency  in  the  lone  pair  density,  in  shape,  size,  position,  and  hydrogen  bonding 
geometry  suggested  that  the  electrostatic  potential  surrounding  the  oxygen  atoms  might  be 
relatively  consistent  in  shape  and  magnitude.  This  however  was  not  the  case.  Fig.  7  demonstrates 
the  electrostatic  potential  of  the  0(1)  hydroxy  group  lone  pair  region.  The  significant  portion  of 
the  negative  ESP  extends  approximately  180°  around  the  oxygen.  A  maximum  ESP  of  ~  -0.20 
eA4  was  observed  for  this  hydroxy  group.  The  additional  negative  regions  displayed  represent 
the  potential  created  by  the  n  density  of  the  aromatic  ring.  The  oxygen  atom  forms  a  continuum 
with  the  negative  potential  above  and  below  the  aromatic  ring. 


Figure  7.  Electrostatic  potential  isosurface  around  the  0(1)  atom.  The  blue  surface  represents 
+1.0  eA'1,  while  the  red  surfaces  represent  -0. 15  eA'1. 

Before  discussing  the  0(2)  hydroxy  group,  it  is  important  to  recognize  the  orientation  of  the 
hydrogen  atom  H(20).  A  hydroxy  group  attached  to  an  aromatic  ring  almost  exclusively  has  die 
hydrogen  atom  in  the  plane  of  the  ring.  Obviously,  a  hydroxy  group  bound  to  an  sp3  carbon  can 
freely  rotate.  In  the  solid  state,  the  orientation  of  the  OH  group  will  almost  always  be  defined  by 
the  hydrogen  bonding  scheme.  Fig.  8  displays  the  D-ring§  of  the  E2  molecule  and  the  torsion 
angle  (y)  of  49°  between  H(20)  and  H(1 7). 


§  The  C(13)-C(14)-C(15)-C(16)-C(17)  ring. 
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Figure  8.  The  0(2)  hydroxy  group  geometry  showing  torsion  angle  (y)  between  the  H(20)  and 
H(  17)  atoms.  The  0(2)  atom  is  shown  in  red. 


The  negative  electrostatic  potential  around  0(2)  has  much  greater  maximum  value  than  that  of 
the  0(1)  hydroxy  group,  reaching  a  value  of  ~  -0.35  eA'1.  In  addition,  the  ESP  does  not  evenly 
cover  die  lone  pair  region  as  on  the  0(1)  hydroxy  group.  A  shift  toward  the  lone  pair  opposite 
H(17)  appears  to  have  taken  place  (Fig.  9).  The  coverage  area  of  the  negative  potential  seems 
smaller,  covering  near  120°,  versus  the  approximate  180°  coverage  on  0(1). 


Figure  9.  Electrostatic  potential  isosurface  around  the  0(2)  atom.  The  blue  surface  represents 
+1.0  eA'1,  while  the  red  surfaces  represent  -0.25  eA'1. 

It  is  believed  that  the  agonistic  or  the  antagonistic  responses  of  different  estrogens  can  be 
related  to  their  electrostatic  potential  [3],  We  intend  to  construct  a  database  with  ESP  calculated 
for  estrogens  with  different  biological  activities.  Once  enough  data  are  available,  we  hope  to 
identify  specific  topological  features  in  these  estrogens  that  distinguish  binding  ability  of  these 
molecules  to  the  estrogen  receptor. 
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Conclusions 

This  is  one  of  the  first  charge  density  studies  accomplished  for  the  big  (56  atoms)  organic 
molecules  in  the  non-centrosymmetric  unit  cell.  Our  study  showed  that  it  is  possible  to  get  a 
quality  data  for  such  kind  of  system;  the  optimal  coordinate  system  and  slow,  initially 
constrained  refinement  are  the  necessary  tools  for  this  type  of  study.  Since  most  estrogens  are 
identical  in  the  core  structure  of  the  molecule,  this  study  will  provide  a  better  starting  point  for 
future  multipole  refinements  of  other  estrogen  derivatives.  Overall  our  results  pertaining  to  the 
urea  molecule  reasonably  agree  with  the  previous  literature  results  showing  the  quality  of  our 
data. 

For  the  E2  molecule,  the  lone  pair  densities  of  the  oxygen  atoms,  which  are  an  extremely  small 
features  of  the  system,  were  found  to  be  quite  similar.  This  indicates  that  the  oxygen  atom  bound 
at  the  C(3)  position  do  not  have  a  strong  interaction  with  the  #  orbitals  of  the  aromatic  ring  in  the 
ground  state  configuration.  Despite  this  consistency,  the  electrostatic  potentials  in  the  region  of 
the  hydroxy  groups  differ  greatly  in  both  shape  and  magnitude.  Future  studies  may  reveal  how 
the  H(17)  -  C(17)  -  0(2)  -  H(20)  torsion  angle  affects  the  electrostatic  potential  of  the  C(17) 
hydroxy  group.  Obviously  many  more  estrogen  systems  must  be  completed  before  the 
comparative  portion  of  the  study  can  take  place  and  real  conclusions  can  be  made  on  the 
relationship  between  the  electrostatic  potential  and  estrogenic  action. 

We  would  like  to  thank  Dr.  Samuel  Brooks  for  his  contribution.  We  thank  Dr.  A. Volkov 
and  Prof.  P.  Coppens  for  making  TOPXD  program  available  for  us.  We  would  also  like  to  thank 
the  Department  of  Defense  for  their  financial  support  (Grant  #  DAMD 17-00- 1-0468). 


Supplemental  Tables  and  figures 

SI-  experimental  protocol 

S2  -  atomic  coordinates 

S3,  S4  -  atomic  thermal  parameters 

55  -  bond  lengths 

56  -  bond  angles 

S7,S8,S9  -  multipole  populations 
S 10,  S 1 1  ,S  12  -  topological  results 

Fig.  Dl.  Local  multipole  coordinate  system  for  the  urea  molecule 
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